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A microscopic model of non-reciprocal optical effects in an- 
tiferromagnets is developed by considering the case of Cr2 03 
where such effects have been observed. These effects are due 
to a direct coupling between light and the antiferromagnetic 
order parameter. This coupling is mediated by the spin-orbit 
interaction and involves an interplay between the breaking of 
inversion symmetry due to the antiferromagnetic order pa- 
rameter and the trigonal field contribution to the ligand field 
at the magnetic ion. We evaluate the matrix elements rele- 
vant for the non-reciprocal second harmonic generation and 
gyrotropic birefringence. 
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I. INTRODUCTION 

The study of optical phenomena in magnetic sub- 
stances has always been an interesting area of research. 
Optical effects exhibited by magnetic substances can be 
classified broadly into two categories, reciprocal and non- 
reciprocal j^]. Reciprocal optical phenomena are those 
that cannot distinguish between magnetic states that are 
related to one another by time reversal. Typical exam- 
ples of such phenomena are those that involve scatter- 
ing of light by magnetic excitations. Non-reciprocal phe- 
nomena, on the other hand, can distinguish between two 
magnetic states that are related to each other by time 
reversal. A classic example of such a phenomenon is the 
Faraday effect discovered by Faraday in the last century. 
The Faraday effect in ferromagnets, for instance, appears 
as a rotation of the plane of polarization of light which is 
incident along the axis of magnetization. This rotation is 
non-reciprocal in the sense that it changes sign on revers- 
ing the direction of magnetization (which is equivalent to 
time reversal) . Such an effect can be understood in terms 
of an interaction between the internal molecular field of 
the ferromagnet and the incident electromagnetic radi- 
ation mediated by the spin-orbit coupling Clearly, 
such an effect would not only be exhibited by ferromag- 
nets but by any substance that has a net (non-zero) mag- 
netic moment such as ferrimagnets or paramagnets where 
a magnetic moment is induced by the application of an 
external field. 

It is then but natural to ask what happens in the case 
of antiferromagnets. At first sight, it would seem that it 
is impossible for light to distinguish between states that 
are related by time reversal as, in this case, there is no 



net magnetic moment (i.e., the total molecular field is 
zero unlike in ferromagnets). However, this is not true. 
In fact, it has been known for a long time from symme- 
try arguments that certain classes of antiferromagnets 
can show a variety of non-reciprocal phenomena 
although till recently, there have been no reports of ex- 
perimental observations of such phenomena. The util- 
ity of such experiments, where possible, can hardly be 
overemphasized since experiments such as Raman scat- 
tering etc., only probe the antiferromagnetic structure 
indirectly by coupling to the magnetic excitations rather 
than the magnetic ordering itself 

In this context therefore, the discovery of non- 
reciprocal optical effects below the Neel temperature, T/v, 
in optical experiments on Cr203 comes as a break- 
through in the study of antiferromagnetic ordering by 
light. Krichevtsov et al. [|j reported the first experi- 
mental observation of spontaneous non-reciprocal rota- 
tion and circular dichroism below the Nccl temperature 
of Cr2 03. Fiebig et al. Q found that antiferromagnetic 
domains could be observed directly by non-reciprocal sec- 
ond harmonic generation (SHG), leading to the first pho- 
tographs ever of antiferromagnetic domains These 
experiments show that light can indeed couple directly to 
the antiferromagnetic order parameter, thereby leading 
to non-reciprocal effects. As mentioned earlier, though 
such a coupling was anticipated from symmetry consid- 
erations, no microscopic mechanism has been presented 
so far. 

In this paper, we present, in detail, a microscopic 
mechanism that describes all non-reciprocal optical ef- 
fects by considering the case of Cr2 03 where such effects 
have been observed experimentally. In an earlier paper 
|l0| , we showed that all non-reciprocal effects in Cr203 
can be explained by the fact that electric dipole transi- 
tions are allowed in Cr203 below the Neel temperature. 
Here we highlight the actual evaluation of the matrix 
elements relevant for the non-reciprocal non-linear sus- 
ceptibilities. The evaluation of the matrix elements is 
performed within a new cluster model for Cr203, which 
contains the full crystal symmetry of Cr203 and which 
allows the orders of magnitude of all matrix elements con- 
tributing to the non- reciprocal phenomena in Cr203 to 
be predicted. We apply the microscopic model to the ob- 
served phenomenon of SHG 0| and explain how antifer- 
romagnetic domains can be distinguished experimentally. 
We also apply our model to another non-reciprocal effect 
seen experimentally in Cr203 viz., gyrotropic birefrin- 
gence Pi and solve the long-standing question regarding 
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its magnitude. Our model can also be used to describe 
non-reciprocal optical effects in other antiferromagnets 
where inversion symmetry is broken below T/v- 

The paper is organized as follows. In section II, we 
illustrate how symmetry arguments can be used to study 
non-reciprocal effects by considering the macroscopic the- 
ory of SHG in Cr203. In section III, we present the mi- 
croscopic theory of SHG which allows us to obtain all 
polarization selection rules and magnitudes of the non- 
linear susceptibilities. In section IV, we use our micro- 
scopic model to explain the phenomenon of gyrotropic 
birefringence and the associated magnetoelectric effect. 
A summary of our results and conclusions are presented 
in section V. 

II. MACROSCOPIC THEORY 

A clue to the origin of non-reciprocal effects can be 
obtained from macroscopic (symmetry) considerations of 
the susceptibility tensors. The optical properties of a 
medium are characterized by linear and nonlinear sus- 
ceptibility tensors. If the susceptibility tensors are known 
for a given medium up to a certain order n, then, at least 
in principle, the nth-order nonlinear optical effects in 
the medium can be predicted from Maxwell's equations. 
Physically, the susceptibility is related to the microscopic 
structure of the medium and can be properly evaluated 
only by doing a full quantum-mechanical calculation as 
we will show in the next section. Nevertheless, one can 
get some information about the susceptibilities just from 
symmetry considerations as dictated by Neumann's prin- 
ciple, which states prj l that any symmetry which is ex- 
hibited by the point group of the crystal is possessed by 
every physical property of the crystal. Then, to investi- 
gate the effect of crystal symmetry on the components 
of the susceptibility tensors, it is necessary to enforce 
the requirement that the tensor be invariant under all 
the permissible symmetry operators appropriate to the 
particular crystal class. Accordingly, some tensor ele- 
ments are zero or are related to others, thereby reducing 
the total number of independent tensor elements. The 
application of Neumann's principle to four-dimensional 
space-time remains valid only for static properties. For 
dynamic processes, like non-reciprocal SHG, where there 
is a preferred direction of time, one considers only those 
symmetry operations that do not include time reversal 
in classifying the allowed tensors, i.e., only the spatial 
symmetry is used to simplify the form of the suscepti- 
bility tensors. Though one can still classify the tensors 
as i and c (time-symmetric/antisymmetric) tensors, one 
cannot use this classification to obtain those that are al- 
lowed. To illustrate these notions, let us now consider 
the macroscopic theory of SHG in Cr2 03 Q. We shall 
consider in detail the case where the spins are oriented 
parallel to the crystallographic z axis and laser light prop- 



agates along the z direction and then we shall generalize 
to the case of arbitrary direction for the orientation of 
the spins and for the propagation of the laser light. 

A. Spins and laser light parallel to z axis 

Above the Neel temperature, Tjv — 307^, Cr203 crys- 
tallizes in the centrosymmetric point group D^d- The 
four Cr'^+ ions in the unit cell occupy equivalent c posi- 
tions along the 3z (optic) axis. Since this structure has 
a center of inversion, parity considerations forbid electric 
dipole transitions in SHG. Magnetic dipole transitions 
that are related to the existence of an axial tensor of odd 
rank, or electric quadrupole transitions related to the ex- 
istence of a polar tensor of even rank are however allowed. 
Below T/v, the spins of the Cr'^+ ions orient along the z- 
axis. The spin ordering breaks time reversal symmetry, 
R, and as SHG is a dynamic process, only symmetry op- 
erations of the crystal that do not include R may be used 
to classify the allowed tensors for the susceptibilities. For 
Cr203, the remaining invariant subgroup is D^j. New ten- 
sors are allowed in this point group, for instance, a polar 
tensor of odd rank, that allow electric dipole transitions 
in SHG. 

From Maxwell's equations, the source term corre- 
sponding to SHG can be derived by considering the 
contributions to (nonlinear) magnetization, M^^'^) = 
eoC 7(2'^) : E(") E(") and polarization, P^^'^) = e^x^^'^' : 
E^'^' 'E'--^\ in the point group D^. In general, one 
should also consider contributions to nonlinear suscep- 
tibilities related to electric quadrupole transitions, viz., 

Q(2c^) ^ ^i£lo^(2t^) . gi^(,g ^j^g l_,j^g-g physics 

associated with such transitions is not very different from 
the magnetic dipole transitions, we will in the following 
concentrate on the electric and magnetic dipole contri- 
butions to the nonlinear susceptibilities. If we assume 
that laser light propagates along the optic axis, there is 
only one independent component 7™ = jyyy — —jyxx — 
~lxyx = ~lxxy of the nonlinear magnetic susceptibility 
(axial tensor of third rank) and analogously for the non- 
linear electric susceptibility (polar tensor of third rank), 
Xe = Xyyy = -Xyxx = -Xxyx = -Xxxy The source term 
S(r,i) in the wave equation, 

[Vx (Vx) + (l/c2)aVa<2]E(r,<) = -S(r,i), (1) 

can be written in a dipole expansion as [ p^ , 

S(r,0 = (2^ + VxiMM 

-v^ + ...) . 

If we now assume that E, P and M can be decomposed 
in a set of plane waves and consider a circular basis with 
E = £'+(3+ + i?_e_ + Eziz, where e± — -y={exiiey), one 



2 



obtains (note the different basis choice with respect to 
refs. § and 0]), 

S = I ^- j = I (7™ +^Xe)El j , (3) 

where huj is the energy of the incoming light beam. Note 
that incoming right circularly polarized light leads 
to left circularly polarized light and vice versa in 

SHG. Above T^, the electric dipole contributions dis- 
appear (xe = 0) and therefore the SHG intensities I± 
are identical while below T^r, (xe 7^ 0) and the intensi- 
ties I± oc 1 7m ± iXel'^E^ are different for right and left 
circularly polarized light, as observed experimentally 0. 
As the tensor x appears only below Tjv, it is natural to 
assume that it is proportional to A(T), the antiferromag- 
netic order parameter. Non-reciprocal SHG can then be 
understood as arising from an interference between 7^ 
and Xe- 

B. Spins and laser light in arbitrary directions 

An extended study of SHG in Cr2 03 can be carried 
out by considering the propagation of light and the ori- 
entation of the spins in the crystal to be arbitrary, the 
motivation of such a study being the fact that application 
of an external magnetic field causes the spins to orient 
along different directions. The spin-flop phase which oc- 
curs in Cr203 when a static magnetic field is applied 
along the z direction can be studied using these results. 
In Table I, we present the results of the macroscopic sym- 
metry analysis below Tjv for light propagating along the 
z, y and x directions and for spins oriented parallel to 
z, y and x axis. The three components of the source 
term are indicated for the various cases. Note that, in 
the first row which corresponds to the case that the spins 
are aligned parallel to the z axis, the interference effect 
occuring for light propagating along the z direction is ab- 
sent when light propagates along the y and x directions, 
as has been verified experimentally [ p^ . 

When the spins are aligned along the z, y and x direc- 
tions the invariant point subgroups of crystal symmetries 
are D3, Cih and C2 respectively. Table I contains the 
results for the source term components for all three pos- 
sibilities. Table H contains the components of the source 
term corresponding to electric quadrupole transitions for 
light propagating along the z, y and x directions and 
spins aligned parallel to the z axis. 

III. MICROSCOPIC THEORY OF SECOND 
HARMONIC GENERATION 

The macroscopic theory discussed so far is based purely 
on symmetry considerations. It can neither provide mag- 



nitudes of the components of the magnetic (7) and elec- 
tric (x) susceptibilities nor can it specify the details of 
the interaction and the origin of the coupling between 
light and the antiferromagnetic order parameter. Both 
these can be obtained only from a microscopic approach 
which is described below. In what follows, we consider 
in detail the case when the spins are aligned parallel to 
the z axis and light propagates along the same direction. 
This is the most representative case for the observation of 
non-reciprocal phenomena in Cr203. Elsewhere, we have 
shown that non-reciprocal effects can be understood by 
the fact that electric dipole transitions are allowed below 
Tjv [0. We use this in calculating the nonlinear sus- 
ceptibilities of Cr203. The method of calculating non- 
linear susceptibilities of dielectric media was developed 
by Armstrong and coworkers In this approach, the 
current density induced in the system by incident elec- 
tromagnetic radiation is calculated using semi classical 
perturbation theory. 

Let us consider the electromagnetic field to be a super- 
position of three harmonic waves whose frequencies coi, 
i = 1,2, 3, satisfy the condition W3 = wi +uj2- We choose 
the vector potential of this field to be 

A(r, t) = diQi exp i(k • r — ujit) + h.c. . 

i 

where Hi is a unit polarization vector and qi determines 
the strength of the electric field, Ei = qiUJi/2 (where 
the Coulomb gauge is being assumed). The perturbation 
is then described by the usual form of the interaction 
hamiltonian 

g 

— Aj • pj + - — Aj • Aj . 

j 

Assuming for simplicity that the ground state wave func- 
tion $ is a product of one-electron wave functions, $ — 
H„(/)„, the induced current density is given in the inter- 
action representation by 

V('/'„|ev|(/)„) = V(0„|p - eA|0„)— . 

^-^ nie 

n n 

Nonlinear susceptibilities are now calculated to the de- 
sired order in field strengths. We are in particular inter- 
ested in terms that are quadratic in fields and the source 
term at = wi -I- W2- Thus we only consider those terms 
in the perturbation expansion for the current density that 
are proportional to (71(72, viz., proportional to E1E2 and 
with time dependence exp(±ict;3t) and write 

(ev(w3)) = (ev(+)(w3)exp(iw3t)) 

-I- (ev(")(tj3) exp(-iw3i)) . 

By multiplying the terms proportional to exp(±iw3t) by 
exp(=Fzk3 • r) and then taking the expectation value. 
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we get the induced current density corresponding to the 
source term at lu^ and at ki + k2 = ka as 

(J(k3,W3)> = (*|ev(-)(cj3)expi(k3T-Cc;30l*) ■ (4) 

One can now obtain the expression for the nonhnear sus- 
ceptibihties by using the relation J = fiodP/dt as also 
the defining relations for the susceptibilities that relate, 
for example, P and E. 

We now turn our attention to the specific case of 
Cr203. As explained in the previous section, we are 
interested in the nonlinear polarization p(2") and the 
magnetization M'^^'^^ The nonlinear susceptibility, x is 
only allowed below the Neel temperature where spatial 
inversion symmetry of the crystal is broken by antifer- 
romagnetic ordering. In an earlier paper [p^ , we have 
shown that this leads to electric dipolc transitions be- 
tween the **A2 and '*T2 levels of Cr203. Thus, below T/v, 
one can use the electric dipole approximation with the 
interaction hamiltonian of the form 



H = 



A simple calculation now yields the result for the induced 
current, 



(J(k3,W3)) = ^ 51 ^ 

|m),|n) 

{<^\J{k3,0J3)\m){m\Ef,rf,\n){n\E,r,\'i>) 

{uJm - 2u;){LJn - Uj) 

^ ($|^^r^|m)(m|J(k3,c^3)|n)(n|g,r,|$) 

{uJrn + W)(w„ - LO) 

{^\Ef,r^\m){m\E,rM{n\3{k3,iU3)m 
(ujn + 2a;)(u;,„ + lu) 



(5) 



where E^^ are the components of the incident electric 
fields and the definition is to be used. In the above 
expression we have also taken (as in SHG) lui — lu2 — lu 
and uj^ = 2oj. All factors exp(ik • r) have been set to 
unity (electric dipole approximation). The problem now 
reduces to that of evaluating the appropriate transition 
matrix elements in the above expression with the Cr^+ 
wave functions. To do this, we consider a cluster model 
that has the correct symmetry of Cr203, viz., D^d- This 
model [ p^ contains only two Cr ions in a magnetic unit 
cell whereas in the actual structure of Cr203 there are 
four. While this would lead to quantitative differences, 
the basic physics remains the same. 

The ground state wave function of the magnetic ion can 
be written as fl^ |<i>) = |t; 
states |e) = le^-^, le^^^) where 



'^\t'^\t'^^) and the excited 
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= 3-i/2[\/2|42_^2) + 14.) 
= -3-i/2[V2|rf,,) 



m \pz) 
m \px) 
m \Pv) 



(6) 



eW = 2,-^'^[\dxy) - V2\dyz)] + m\Py) 

e(2) = 3-i/2[|d,2_^2) + V2|4.)] + ^73 1?.) , 

Here the \d) and \p) states are the Cr-3ci and Cr-4p or- 
bitals respectively. In deriving the above wave functions, 
we have chosen the axis of quantization to be parallel to 
the crystallographic z direction. The following conven- 
tion for the d orbitals is chosen. 

dxy = -^^^ I TZ ) 



An 
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Here the R{r) are the appropriate radial functions for 3d 
orbitals. Spin quantum numbers in (^ have been sup- 
pressed for clarity. We also note that the expression (||) 
only includes the effect of the hemihedral part of the trig- 
onal distortion which is the most dominant interaction. 
This interaction is of the form rjz and leads to a mixing of 
Cr-3c? states and Cr-4p states {px,Py,Pz) with coefficients 
771, ri2 and 773 being proportional to 77, the trigonal field. 
In addition to this mixing, the spin-orbit coupling also 
causes a mixing of the ^2 and e orbitals. As explained in 
our previous paper, it is this interplay between the trig- 
onal distortion and the spin-orbit coupling that leads to 
electric dipole transitions below T/y. Consequently, we 
focus our attention on those transition matrix elements 
in (^ that are proportional to both the spin orbit inter- 
action and the trigonal field. It is easy to see that such 
matrix elements are of the form 

{t2\L,S,\e){e\r\m)(m\E^r^\n){n\E,r,\t2) , (7) 

where we have used (7Ti|p|77) — —ime{oJn — ^m){'ni\v\n) , 
etc. Let us now consider the case when the incident light 
is right circularly polarized 0, viz., E = Ee+ where e± = 
-^(ca; ± iey). The relevant transition matrix element is 
now given by 



{tf\L,Sz\e^'^){e^'^\r\e^''^){e^'''^\r+\n){n\r+\tf) 



where the states \n) are higher-energy orbitals of odd 
symmetry. Let us now consider that part of the tran- 
sition matrix element given above which corresponds to 
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the emission of the second harmonic, viz., 
(4^'|L,5,|e(^"))(eW|r|e(^)) . 

Since the excited state Cr orbital |e'--'^) is a mixture of 
the Cr-3ci and 4p states, the emission matrix element has 
contributions of the form 



(8) 



where rj is the appropriate admixture of the p the d or- 
bitals determined by . Eq. || describes the contribu- 
tion to the ED matrix element of a single Cr ion. The 
response of the crystal is given by the coherent summa- 
tion over the contributions of each ion. In Cr203 there 
are two pairs of sites in the unit cell, the A/B and the 
A' / B' sites related by inversion symmetry. Since they 
give identical contributions we need to consider only one 
of the pairs. The contribution of a given unit cell to the 
ED matrix element is then proportional to 



A {tia{Sz)a + r]B{Sz)B) 



(9) 



where A is the spin-orbit coupling constant (see ref. |10[ ) 
and TjA and rjB are proportional to the respective local 
trigonal fields. Inversion symmetry demands that —riB = 
TjA = rj. The total contribution to the ED matrix element 
is then proportional to 

M{Sz)a - {Sz)b) const. = XrjA{T) 

all unit cells 

where A(T) is the antiferromagnetic order parameter. 
One can use the explicit form of the orbitals given by (|^) 
and evaluate the desired matrix elements. For example, 
the contribution to SHG by exciting the electron in the 

(2) 

orbital is proportional to 

\r]A{T)[{di^^2_y2)\Lz\dx:y){dxy\x\py) x 

{Py\z\dyz){dyzm-r)^\dz.)] . 

A similar expression can be written for the contribution 
arising from t^"^ . From our results we find that (a) the 
t'^^ orbital does not contribute to SHG; This is not sur- 
prising as our mechanism of SHG is mediated through the 
diagonal part of the spin-orbit interaction and Lzt'^^ = 0. 
(b) the electric field of the emitted radiation at frequency 
2lu is along the x — iy direction, i.e., incoming right circu- 
larly polarized light generates the second harmonic with 
opposite polarization. Likewise, it can be verified that 
left circularly polarized light generates a second harmonic 
that is right circularly polarized. Thus our microscopic 
results reproduce the selection rules that are given by 
macroscopic theory. This is because the expression (^) 
for the Cr'^'*" orbitals has been obtained by incorporating 
the full crystal symmetry of Cr203. 

The constructive interference in emission of the contri- 
butions to the ED matrix elements below T/v, as given 



by Eq. ||, is a consequence of the symmetry of the spin- 
ordering in Cr203. We reemphasize that this is because 
our mechanism is a one-ion mechanism that invokes the 
diagonal part of the spin-orbit interaction. Photons ab- 
sorbed by a certain ion are, in this mechanism, reemitted 
by the same ion. The phase of the emitted electromag- 
netic wave then depends on the relative direction of the 
spin-ordering with respect to the direction of the trigonal 
distortion at the given site (||). The coherent interfer- 
ence of the waves emitted by all ions may then interfere 
constructively. The effect is non-reciprocal since the di- 
rection of the local trigonal distortion is invariant under 
time-reversal, whereas the spin-order is not. 

Let us now consider the case T > T^. Now, inversion 
is a good symmetry of the crystal and the electric dipole 
transitions vanish identically. In this case, one has to 
consider the non vanishing term of the lowest order in the 
multipole expansion of exp(— iks • r). This is, of course, 
the magnetic dipole term which can be written as 



'^'"e"- \Tn).\n) 

(m| cxp{-iki-r)aif_,Pf_,\n){n\ exp(-zk2 •r)a2,yp,y|$) 
exp(iLJ3<) 

i _L 

(uJrn - L^n + ^i) (w„ + CJ2) 



(10) 



The calculation of transition matrix elements in this 
case is much simpler as now the dominant contributions 
contain neither the spin-orbit interaction nor the trigo- 
nal distortion. It is easy to verify that one only has to 
consider contributions to the transition matrix elements 
that are of the form 



^ ^(i«|L|e(^-))(e(^-)|r+|n)(n|r+|t«) 



(11) 



i=l,2,3i=l,2 



at any given Cr site. (Here again we have assumed that 
the incoming light is right circularly polarized.) The con- 
tributions from different sites can be added up coherently. 
As in the previous case of the ED transition, we have also 
verified using (H) that the correct polarization selection 
rules are reproduced by the above expression below and 
above Tpf. 

The microscopic expression for the dynamical current 
operator together with the Cr ion wave functions (^ en- 
ables us to estimate the magnitudes of the nonlinear sus- 
ceptibilities. Since we do not have the correct wave func- 
tions of the ea;cite(i states, our estimates are approximate, 
correct only to an order of magnitude. In particular one 
may also include contributions from the 0-2p orbitals in 
(^). Such a contribution would not alter the mechanism 
for the ED-transition explained above but would lead to 
certain quantitative corrections. From expressions (|^) 
and (^, as also the definition of the macroscopic suscep- 
tibility X obtain the order of magnitude of x^^^^ 
approximately 
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(12) 



Here, Xq ~ 5000A is the wavelength of the emitted hght, 
flo « 0.69 A is the radius of Cr'^+j Uq is the density of 
Cr ions in CrsOg (~ 3.3 x 10^^ m^^), A « 100 cm~^ is 
the spin-orbit interaction 0,|l8|, E^-Et^ w 8000cm~i is 
the difference in energy between the ^2 and the e orbitals, 
r] « 350 cm-i is the trigonal field Ep - Ed ~ & x 
lO"' cm^^ is the difference in energy between the d and 
the p orbitals that are mixed by the trigonal distortion 
and A(T) is the antiferromagnetic order parameter. The 
temperature dependence of x is solely through that of 
the antiferromagnetic order parameter and consequently 
X vanishes above T^. At temperatures much lower than 
Tat we find that ^ 8.7 x IQ-ie Coul N"!. 

Similarly, we find that the magnitude of ^^^^"f is ap- 
proximately 
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(2c^)| 



\2hujJ Amc 



(13) 



As explained earlier, this is the MD contribution to SHG 
which depends neither on the trigonal distortion nor the 
spin-orbit interaction. It is also independent of temper- 
ature and we estimate it to be of the order of 11 x 10^^^ 
Coul N~^. The contributions to J from the nonlinear 
susceptibilities x and 7 can produce interference effects 
in SHG only if the first two terms in the r.h.s. of equation 
(^ are of the same order of magnitude, wz., \x\ — 
From our estimates of x and 7 above we find that this is 
indeed satisfied thereby leading to an interference effect 
in SHG which is non reciprocal and which vanishes above 
Tn- 



IV. GYROTROPIC BIREFRINGENCE AND THE 
OPTICAL MAGNETOELECTRIC EFFECT 

The ED transition in the optical region that is allowed 
below T/v in Cr203 can also be seen in one-photon ex- 
periments. In this section, we consider the phenomena of 
Gyrotropic Birefringence (GB) and the associated Opti- 
cal Magnetoelectric Effect that are one-photon processes 
in contrast to SHG which is a two-photon process. 

GB is a non-reciprocal optical effect that appears as 
a shift in the principal optical axis along with a change 
in the velocity of propagation of light. The possibility of 
observing this effect in Cr203 and using it to distinguish 
between antiferromagnetic domains was first pointed out 
by Brown and co-workers jj] . Their analysis of the prob- 
lem was purely macroscopic and they showed that GB 
appears in the form of a polar c tensor (one that changes 
sign under both space and time inversion) below T^v. 



Later, Hornreich and Shtrikman |l9| presented the first 
quantum mechanical treatment of this problem. They 
showed that the gyrotropic birefringence tensor can be 
described in terms of electric quadrupole and magneto- 
electric effects. From their results, they estimated that at 
optical frequencies, the induced rotation of the principal 
optical axes of Cr203 would be of the order of 10~® — 10~^ 
rad and the magnetoelectric susceptibility would be of 
the order of 10^® (in dimensionless units). This phe- 
nomenon was also analyzed by Graham and Raab pc| ] 
using a multipole theory of wave propagation. 

Recently however, Krichevtsov et al. Q| observed for 
the first time, spontaneous non-reciprocal rotation of the 
optical axes in Cr203. They found that the observed ro- 
tation and the magnetoelectric susceptibility were 4 or- 
ders of magnitude larger than those predicted by Horn- 
reich and Shtrikman. They also found that the observed 
temperature dependence of these non reciprocal effects 
corresponds roughly to that of the antiferromagnetic or- 
der parameter, something that does not follow obviously 
from previous calculations [|9| . The intensity of these ef- 
fects and the temperature dependence led them to specu- 
late that these effects were attributable to ED transitions 
in the optical range. We now show that this is indeed the 
case and the mechanism we discussed earlier does lead to 
effects that are of the same order of magnitude as those 
observed. 

To see this, we calculate the magnetoelectric suscep- 
tibility in the optical region, allowing for an ED transi- 
tion below Ttv. The calculation is in many ways simi- 
lar to that presented earlier for the nonlinear suscepti- 
bilities. The magnetoelectric susceptibility is defined as 
M^'^' = a : E*^") and we estimate it by calculating di- 
rectly, the quantity (M) = g^B(L + 2S) in perturbation 
theory. Since in this case we are only interested in tem- 
peratures below the Neel temperature, we use the electric 
dipole approximation and consider the linear response of 
the system. It is easy to verify that the induced magne- 
tization is given by 



(M) = jRe[V($|M|m)^^^^^-^^exp(-zc^t) 

\m) 



($|M|m)^^ ^ exp +zc^t 

(CJ +UJrnj 



where we have assumed the electric field to be of the form 
E = Eo cos(k-r— wi). Choosing (as in the experiment |^) 
the incident light along the x direction and the electric 
field E|ji , it is easy to see that a typical matrix element 
that would contribute to the magnetoelectric susceptibil- 
ity is of the form, 

{t2\L,\e){e\L,S,\t2){t2\zE,\t2) . 

Note that {t2\zEz\t2) would be proportional to the trig- 
onal mixing of the 3d and the 4p orbitals. A straightfor- 
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ward calculation using to evaluate appropriate transi- 
tion matrix elements gives us the expression for the mag- 
netoelectric susceptibility, 



VI. APPENDIX 



A 



4/ioCe 



h{uj - LUm) 



Ee — Ep — Ed 



2o A(T) 



(14) 



in dimensionless units. Here, Uq is the density of Cr ions 
in CraOa (~ 3.3 x lO^^ m'^) and h{Lu - w„) ~ 0.5 eV 
in the region of experimental interest. Thus by evaluat- 
ing the above expression, we estimate axx ~ 0.2 x 10~"* 
which is of the same order of magnitude as that observed 
experimentally. This also means that the non-reciprocal 
rotation would be ~ 10~* rad. Since the ED process 
that we consider couples light to the order parameter, 
the observed temperature dependence follows naturally 
from our theory. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have developed a microscopic the- 
ory of non-reciprocal optical effects observed below Tn 
in Cr2 03. We have shown that these effects can be ex- 
plained by an electric dipole process that arises from an 
interplay between the spin-orbit coupling and the trigo- 
nal distortion of the ligand field. Such a process couples 
light directly to the antiferromagnetic order parameter. 
In contrast to other ED mechanisms that have been con- 
sidered in the literature so far [ p6| , our mechanism is 
a one-ion mechanism that couples light to the sublat- 
tice magnetization rather than any magnetic excitations. 
Photons absorbed by a certain Cr'^+ ion are also emitted 
by the same ion. The coupling to the antiferromagnetic 
order parameter then occurs throught constructive inter- 
ference of the photons emitted by different Cr-ions. 

It might seem that our mechanism leads to effects that 
are weak as they involve the trigonal distortion and spin- 
orbit interaction, both weak effects by themselves. How- 
ever their interplay leads to an electric dipole transition 
whose oscillator strength is large in the optical region. 
Consequently, observed effects are strong. Though we 
have used this mechanism to explain successfully the 
phenomena of Second Harmonic Generation, Gyrotropic 
Birefringence and the Optical Magnetoelectric effect that 
have been observed experimentally in Cr203, our theory 
can be generalized to all materials where i) the magnetic 
ion is not at a center of inversion and ii) inversion is still 
a macroscopic symmetry above but is broken below 
Tjv due to the ordering of the magnetic ions. Thus, non- 
reciprocal effects should be observable, for example, in 
the cuprate Gd2Cu04 below the ordering temperature of 
the Gadolinium magnetic subsystem, TAr(Gd)=6.5K pl| ] 
when inversion symmetry is broken as also in V2O3 and 
MnTiOa. 



The components for the nonlinear susceptibilities (x) 
and (7) in Table I are as follows: 



Xxxy 



Xe 


= Xyyy = 


Xyxx — Xa 


Xei 


= Xyyy 




Xe2 


= Xyxx 


= Xxyx 


Xa 


= Xyxz — 


~Xxyz 


Xai 


= Xyxz 




Xa2 


= Xxyz 




Xb 


= Xyzx — 


~Xxzy 


Xbi 


= Xyzx 




Xb2 


= ~Xxzy 




Xi 


= Xzzy = 


Xzyz = Xyzz 


X2 


= Xzyx 




X3 


— Xxxx 




X4 


= Xxyy — 


Xyyx — Xyxy 


X5 


= Xzxx — 


Xxxz — Xxzx 


X6 


= Xzyy = 


Xyzy — Xyyz 


X7 


= Xzzx — 


Xzxz — Xxzz 


Xs 


= Xzzz 




X9 


= Xyyz — 


Xyzy — Xzyy 


Xio 


= Xzxy 




7m 


= lyyy ~ 


^yxx — 7a^; 


7mi 


= lyyy 




7m2 


— ^xxy 


— ^xyx — ^ 


71 


= Izzy = 


Izyz — lyzz 



^xxy 



Above the Neel temperature all x's and 71 vanish (and 
7mi = 7m2 = 7m)- ^ is then reasonable to assume that 
below Tjv 7i is much smaller than all other matrix ele- 
ments since it is of magnetic dipole character and only 
allowed due to the breaking of inversion symmetry below 
the Neel temperature. The components of the electric 
quadrupole susceptibilities in Table II are as follows. 

Xl = Xxxyy — Xyyxx 
X2 = Xxxzz — Xyyzz 
X3 = Xzzxx = Xzzyy 
X4 = Xzxyy — Xzxxx — Xzyxy 

— ~Xzyyx — ~Xxyzy = ~Xxyyz — ~Xxzyy 
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direction of incident light 




z 


y 


X 


z 


( -f^{El-El)-2xeE^Ey \ 
1 -Xe{El - El) - 2-fmE.Ey 1 


I ' \ 

1 -(Xa + Xb)E:,E, - XeEa^Ex j 


( (Xa + Xb)EyE, \ 

\ ImEl 1 


y 


/ (X3 + lm,)El + (X4 - lm,)El \ 

HXi - lm2)ExEy 

\ X5EI + x^El ) 


/ XsEl + X7EI + 2X5E^E, \ 
\ X^El + xsEl + 2x7E,E, J 


( XiEl + xiEl \ 

2(X9 - li)EyE, 

\ (X9 + lrm)El + (X8 + ll)El j 


X 


1 Jm,El - ^ra.El - 2Xe,ExEy \ 

-Xe,E^ + Xe,El-2jm,E,Ey 
\ (XIO + X2)E:,Ey J 


^ -Xe^El + xiEl + (Xai + Xb,)E^E, j 


( -{Xa, + Xb,)EyE, \ 

Xe,El+xiEl-2^^EyE, 
\^m,El + ^iEl + 2xiEyEj 



TABLE I. Predictions for the componentsof the source term S from macroscopic symmetry considerations. A constant 
term Alu jc has here been omitted for simplicity. The three cohimns are the predictions for light incident in z, y and x 
directions respectively. The three rows are for the moment of the Cr^""" spins ahgned parallel to the crystallographic z, y 
and X axis respectively. 7i (i = m, mi, m2, 1) and Xi (* = e, ei, 62, o, oi, 02, 6, 61, 62, 1, 2, 3, 10) correspond to the remaining 
independent components for the magnetic and electric susceptibilities respectively. See the appendix for an explicit account of 
these components. 





direction of incident light 


z 


y 


X 




(UEl-El)\ 




( X2EI + XiEl \ 


z 


-2x4ExEy 


(^XiE'x + X2E'.^ 


-2x4EyE, 




\UEl + El)j 




V -X.EI J 



TABLE IL Electric quadrupole contributions to the source term, S. The constant AuP' !<? has been omitted for simplicity. 

The three columns are the predictions for light incident in 2:, y and x directions respectively and for the Cr^"*" spins aligned 
parallel to the z axis. Xi (* = !> 2, 3, 4) correspond to the remaining independent components. (See the appendix). 
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